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SUMMARY 

Experimental heat -transfer and drag coefficients were obtained for 
ethanol drops burning in the presence of low-frequency chamber-pressure 
fluctuations of predominant frequencies of 90 to 100 cps. These coeffi- 
cients were within ±10 percent of values determined from transfer rela- 
tions obtained previously for relatively steady combustion in rocket 
combustors. Also, the theoretically predicted effect of completeness 
of fuel-spray vaporization on combustor performance was compared with 
experimental results. The percent by volume of the total fuel vaporized 
was within ±2 percent of the square of the characteristic-exhaust- 
velocity efficiency of the combustor. 


INTRODUCTION 

Knowledge of the rates of exchange of heat, mass, and momentum 
within combustion chambers is needed to determine relations existing 
between the acceleration of combustion-gas streams and the vaporization 
rates of fuel sprays accelerating and burning in the combustor . In 
reference 1, heat-transfer and drag coefficients were determined for 
ethanol drops burning in a relatively stable combustion environment. 
Such information has proved useful for predicting the combustion per- 
formance of rockets (ref. 2). The combustion of liquid propellants in 
rocket engines, however, is often accompanied by pressure fluctuations 
in the combustion chamber as well as in the propellant injection system. 
Similar heat -transfer- and drag-coefficient data for drops burning in 
the presence of combustion oscillations are therefore needed to obtain 
a more general understanding of the combustion process. 

This investigation was made to determine the effect of an environ- 
ment of low-frequency combustion oscillations on the vaporization of 
ethanol drops burning in a rocket chamber. Drop-size distribution and 
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drop-velocity data were obtained for ethanol sprays burning in the pres- 
ence of low-frequency fluctuations of gas pressure by photographing the 
spray in a low-thrust rocket combustor with a high-speed camera developed 
at the Lewis Center (ref. 3). Thus, from these data and from combustion- 
gas-velocity data, it was possible to calculate heat-transfer and drag 
coefficients, which were compared with the values obtained for steady- 
state combustion (ref. l). 


RESULTS AND DISCUSSION 
Experimental Investigation 

The rocket combustor, high-speed camera, and auxiliary equipment 
that were used in this investigation are shown in figure 1. A descrip- 
tion of the development of the photographic technique and the procedure 
for operating the combustor are given in reference 3. With this equip- 
ment, it was possible to obtain size-distribution data for drops burning 
in a combustion chamber. 

In order to vary the fineness of atomization, an injector was used, 
as shown in detail in figure 2. Ethanol (95 percent ethanol and 5 per- 
cent water) was injected into the rocket chamber through the center 
tube. A pair of nitrogen-gas jets was impinged upon the fuel jet at an 
angle of 45° and at a distance of 1 inch from the fuel orifice in the 
injector face. Using a nitrogen-gas flow rate of 0.027 pound per second 
produced a relatively fine spray. For the second set of test conditions, 
a coarse spray was produced with a nitrogen-gas flow rate of only 0.002 
pound per second. Liquid oxygen was injected into the rocket chamber 
through the two parallel rows of orifices in the face of the injector 
shown in figure 2. Also, a nitrogen-gas flow rate of 0.07 pound per 
second was used for internal cooling of the combustor windows shown in 
figure 1. 

Traces of chamber -pres sure variation with time during the rocket- 
firing test (fig. 3) were obtained with a strain-gage pressure transducer 
and a direct-writing oscillograph for two atomization conditions. Rela- 
tively low-frequency pressure fluctuations were observed in both cases. 
The amplitudes of the oscillations as well as the time-averaged chamber 
pressure, however, were considerably greater for combustion of the finer 
spray. By means of a tape recorder and frequency analyzer, the data 
shown in figure 4 were obtained to determine the frequency spectrum of 
the chamber-pressure fluctuations. The frequencies ranged from 50 to 
200 cps and peaked at approximately 100 cps for the fine spray and ap- 
proximately 90 cps for the coarse spray. Values of APp/P^ were ap- 
proximately 37 percent for the fine spray and 10 percent for the coarse 
spray (fig. 3). (All symbols are defined in appendix A.) 
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The chamber-pressure fluctuations were considerably below the lon- 
gitudinal acoustic frequency of the combustion chamber , which was ap- 
proximately 1200 cps. Also, they could not be attributed to combustion- 
gas residence time (approximately 2 msec). The fluctuations of the 
chamber pressure, however, were reproducible and, therefore, useful to 
this investigation even though their source was not determined. 

Photomicrographs of the fuel spray were obtained at camera stations 
6, 10, and 15 inches downstream from the injector. The photomicrograph 
shown in figure 5 was taken at the 10-inch station. Droplet images mag- 
nified 14 times actual size were obtained with the aforementioned camera 
lens system. Drop diameters were measured directly from contact prints 
with an eyepiece having small scale subdivisions of 0.10 millimeter and 
thus were accurate to within ±7 microns. Only drop images that were in 
sharp focus were measured. Table I shows the drop-size-distribution 
data obtained for combustion of relatively fine and coarse sprays at 
distances of 6, 10, and 15 inches downstream from the injector face. Ap- 
proximately 40 photomicrographs were used to compute each set of drop- 
size data. 

Drop-velocity data were obtained with an accuracy of ±15 feet per 
second with a high-speed rotating mirror incorporated in the camera unit 

to track the burning drops and to stop their images on the 9~inch-wide 

infrared film. Since large drops moved at relatively low velocities as 
compared with the high velocity of small drops, a low mirror speed was 
required for them. Several pictures over a range of mirror speeds were 
thus needed to obtain the average velocity for a given-size drop at a 
particular downstream distance. The drop-velocity data (fig. 6) are 
average values obtained from at least 40 photomicrographs. Only one 
photomicrograph was obtained for each test firing of the rocket. 

Average values of combustion-gas velocity were obtained from four 
runs at each station (fig. 7) with a streak -type 16 -millimeter camera 
(100-ft/sec film speed). This camera was used to photograph the flame 
through a combustor window with a l/8- by 1-inch slit located in the 
same plane as, but at an angle of 90° to, the windows used for the high- 
speed tracking camera. Random variations of ±10 percent in the 
combustion-gas velocity were observed for conditions of high -amplitude 
chamber-pressure fluctuations. 


Heat -Transfer Relations 

The rate of heat transfer dQ/dt to a single drop of fuel vaporiz- 
ing and burning in a high-velocity and high-temperature combustion-gas 
stream may be expressed as 


§ - % i = - Dk g " 


(i) 


where is the total heat gained by an incremental weight of liquid 

dm. In a three-step process prior to actual combustion, the incremental 
weight dm gains sensible heat when the liquid is heated from the in- 
jection temperature T^ to the boiling point T g . The liquid is then 
vaporized by gaining the required heat of vaporization H v ; finally, the 
vapor gains sensible heat in being heated to the temperature Ty. Thus, 

% = H v + =p,!< T s - T !> + c p ( v( T y - T s> (2) 

and the Nusselt number Nu in equation (l) is thus evaluated at the 
temperature Ty. 

Equation (l) may be rewritten in terms of the rate of droplet-area 
change as follows: 


dA 

dt 


:rtANu 


( 3 ) 


where A = 4k g AT/p^Hip. In order to apply equation (3) to a spray of 
drops, it may be rewritten as 


dAip 

dt“ 



D=n 


m 


nXNu ^ n 
D=0 


(4) 


Since 


£ riD 2 = N t d1 q 
D=0 


and 


n - Kfp 

D=0 

equation (4) becomes 




dt 


( 5 ) 



The total change in spray area is thus a function of the area change of 
the mean-drop size and a l so the change in the total number of 

drops Nrp due to the complete vaporization of the small drops in the 
spray. 


Since 
ten as 


Nu = Nu 2Q (±5 percent) (ref. 


1), equation (5) may be writ- 


Uu 2 o 


V Z0 

A 



d(ln %) 

d ( ln D lo) 


(6) 


where V^q is the velocity of a drop in the spray having the same di- 
ameter as the mean drop size D 20 and the Nusselt number is defined as 
Nu 20 = W^oA-g* 

The method of evaluating d(ln N T )/d(ln D 20 ) is given in appendix B 
By substituting the value of this ratio into equation (G) and also the 
values of d^D^/dx obtained from the plot of d| q as a function of 
x (fig, 8), it was possible to calculate values of the Nusselt number 


In reference 4 the following expression was derived for the Nusselt 
number of single drops : 


Nu = 2 + 2,58x10 s 


Re Sc gZ 

tt2 


vO.l 



(?) 


In order to obtain a similar correlation for sprays, the following tem- 
perature conditions were used to evaluate Nu^q anc ^ ^ e 20* 

(1) At the liquid-film interface for the fine and the coarse sprays 

at T s of 118° and 111° C, respectively, for both Nu 2Q 

and Re 2 Q 

(2) At the gas-film interface for the fine and the coarse sprays 

at T of 1381° and 1708° C, respectively, for both Nu 20 
S 

and Re 20 

(3) At the average film temperature for the fine and the coarse 

sprays at T f of 750° and 910° C, respectively, for Nu 2Q ; 
at the stoichiometric flame temperature for both the fine and 
the coarse sprays at 5000° F for Re 2 Q (This procedure was 
also used in ref. 1.) 
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(4) At the average film temperature for the fine and the coarse 
sprays at T f of 750° and 910° C, respectively for both 
Nu 20 and Re 20 

Film properties used in the fourth method of computation are given in 
table II, and properties of ethanol and the combustion gases are summa- 
rized in table III. The combustion temperature and gas properties were 
determined by a general method for automatic computation based on equi- 
librium composition of the combustion gases (ref. 5). 

The results of these calculations are shown in figure 9, and it is 
evident that methods 1 and 2 do not agree with single-drop data 
(eq. (7)), whereas methods 3 and 4 do agree with equation (7) to within 
±10 percent. Evaluation of both the Nusselt and the Reynolds numbers 
at the average film temperature (method 4) may be preferable since the 
combustion-gas temperature is treated as a function of the oxidant -fuel 
latio. Vaporization tests made at near-stoichiometric mixtures in a 
combustor would aid in making a choice between these two methods. Since 
the data agree well with equation (7), the chamber -pres sure fluctuations 
appeared to have little effect on the heat-transfer Nusselt numbers for 
the conditions of this investigation. 


Momentum-Transfer Relations 

According to Newton's second law of motion, the force F on a 
vaporizing drop may be equated to the change of momentum of the drop 
with time to give the following expression for the case of varying mass 
and velocity: 


d(mV ) 
dt 


dV dm 

+v 5t 


( 8 ) 


Equating the instantaneous force F on the drop to the aerodynamic 

pressure force of the combustion-gas stream yields the following expres 
sion; 


d(mV *) 1 . .. 

dt ~ 2 p g A d^ V g " C D (9) 

where Cp is the droplet drag coefficient. 

For the case of a spray of drops y equation (9) may be rewritten as 


Z r - D 'V 


D=D, 


D=0 


dt 


m 

-t£ Z 


( 10 ) 


D=0 
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since AV = V g - Vj and A d = itD 2 /4. In order to simplify equa- 
tion (10), the following expressions were derived: 


V 


30 


D=0 


?= D m 


£ 

D=0 


nlT 


and 


( c d)2o( AV )Io 


D=D, 


y nD 2 (AV) 2 C D 


D=0 


W* 2 

L nD 

D=0 


From the experimental drop-size (table i) and the drop-velocity (fig. 6) 
data and the drag-coefficient expressions (ref. 1), these expressions 
were evaluated and found to be accurate to within ±5 percent. Equa- 
tion (10) may be thus simplified as follows: 


d ( V 30 tI T D lo) 

dt 


| %Di 0 (AV) 2 0 (C D ) 20 


(ID 


since 


and 


f 

D=0 


n D° = NmD 


T^30 


f 

D=0 


nD = 


% D lo 


Taking the derivative of the left side of equation (ll) according to 
equation (8) gives 


8 


(VsoNtDIo) _ 3 dV 30 Tr d( NT D3 0 ) 

dt = ^30 — + V 30 at 


( 12 ) 


Equation (12) may be rewritten as 


d ( V 30 W T D 3o) „ „ “30 


dD^ 


dt 


- v 30 n t 


dt 


d(ln N ) d(ln V ) 

II + — 1— + 


d ( ln D 30) d ( ln D 3 o)J 


(13) 


Substituting it into equation (ll) gives the following expression for 
the drag coefficient: 


( c d) 


20 


4 ^ 

3 


V: 


30 


fdD 


30 


( AV )| 0 


dx 


D 


-2 

20 


d(ln H t ) d(ln V 3Q ) 


d(m d|o) 


a ( ln d 3 o | 


(14) 


where V™ is the velocity of a drop in a spray with D = D 


30 

< AV >Z0 = T g - V 20 


30 


and 


and V 


20 


is the velocity of a drop in a spray with 


D = D 


20 ’ 


Values of Nm at each x-distance were determined by the 


method described in appendix B, and corresponding values of D 
obtained from the drop-size-distribution data 


Values of V^q and V 3 q 


were obtained from figure 6 for drops with diameters D = D^q 
D = D ‘ " 


30 
and 


were 


2QJ respectively. 


The relation between V^q and D^q was determined from the slope 
of the curves in figure 10 and the relation between Nrp and D 3 q from 
the slope of the curve in figure 11. By substituting these values and 
the slope dB^/dx obtained from figure 12 into equation (14), it was 
possible to calculate drag coefficients for the mean-drop size D^q. 

In reference 1, the following expression was derived: 


/ -,\-0 2 

fn \ n * -0.84/gA 

( C d) 2 0 “ 0,13 R 20 (=2 


(15) 


Therefore, the drag coefficients determined from equation (14) for two 
film-temperature conditions were plotted as functions of the product 

Re^Q* ®^(gZ/c“ 8 ) ^ (fig- 13). Drag-coefficient data obtained in refer- 
ence 1 are also included for comparison, and agreement of ±10 percent 
with equation (15) is shown. Thus, the presence of chamber -pressure 
fluctuations did not show any appreciable effect on drag coefficients 
obtained in this investigation (fig. 13). 
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Combustor-Performance Calculations 

To determine the performance of the combustor , the characteristic 
exhaust velocity c* was calculated from the expression 


c* = P c g 



(16) 


where A^_ is the area of the nozzle throat, 0.857 square inch, and 

Wm = + Wkt = 0.395 and 0.370 pound per second for chamber pressures 

-L p 1M 2 

of 57.0 and 48.5 pounds per square inch absolute, respectively. 
Nitrogen-gas flow rates were 0.097 and 0.072 pound per second for the 
fine and the coarse sprays, respectively. Results from this study are 
given in table IV in which values of c £heor vere calculated by the 

method given in reference 4. 


The volume percent of the fuel vaporized E was determined from 
the expression 


E 




(17) 


where the product N^D^q was evaluated at each downstream distance x 

by referring to figures 11 and 12. Values of E are plotted as func- 
tions of the downstream distance x (fig. 14), and the total volume per- 
cent of fuel vaporized for each of the two test conditions is given in 
table IV. 


On a theoretical basis, the percent by volume of fuel vaporized may 
be assumed approximately equal to the combustion efficiency, which is 
proportional to the absolute temperature of the combustion gases T g . 

Theoretical rocket -performance calculations show that c */ c £h eor ~ 

Thus, theory predicts that 

\ theor 



Experimental results given in table IV show an agreement with theory of 
±2 percent. 
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SUMMARY OF RESULTS 

Fuel-spray vaporization-rate equations, derived for relatively 
steady-state combustion, may be applied in the case of low-frequency 
chamber-pressure fluctuations (for random amplitude variations as high 
as 37 percent of the average chamber pressure) without modifying the 
expressions; that is, the heat-transfer and drag coefficients based on 
the average chamber pressure and the area-mean drop diameter were within 
±10 percent of the values calculated for steady-state combustion. Also, 
analysis of the combustion data showed that the completeness of vaporiza- 
tion of the fuel spray in the combustor was equal to the square of the 
characteristic-exhaust-velocity efficiency, with an agreement of ±2 per- 
cent. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, July 6, 1962 
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APPEKDIX A 




"D 


P 

c* 

c* 

theor 

D 

D 20 


^30 

dm/dt 

dQ/dt 

E 

e 

P 

f 

g 


SYMBOLS 

area, sq cm or sq ft 

molecular mass diffusivity evaluated at g/(cm)(sec) 

drag coefficient 

root -mean- square velocity of gas molecules, cm/sec 
specific heat at constant pressure, g-cal/(g) (°C) 
characteristic exhaust velocity, ft /sec 
theoretical characteristic exhaust velocity, ft/sec 
drop diameter, microns 

area-mean diameter defined "by general expression: 

(D e ^ ) e “^ = ^2 n ^ >e / X* gives 


D 20 = (E 1 ^ 2 / X n ) 1//2 > microns 


volume-mean drop diameter , f y nD'/ , microns 

vaporization rate, g/sec 

heat-transfer rate 

volume of fuel vaporized, percent 

mean-diameter notation 

force on a vaporizing drop, g-cm/sec^ 

mean-diameter notation 

/ 2 

acceleration due to gravity, 980 cm/sec 

total heat required per unit mass to increase temperature from 
T z to T ' g-cal/g 


H, 


T 
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latent heat of vaporization, g-cal/g 
h heat-transfer coefficient, g-cal/( sec) (cm 2 ) (°C) 

k thermal conductivity, (g-cal) (cm) /(sec) (cm 2 ) (°C) 

l mean free path of air molecules, cm 

M molecular weight 

m mass of liquid, g 

N,p total number of drops in a spray 

Nu heat -transfer Nusselt number for a single drop, hD/k^. 

Nu average heat-transfer Nusselt number, hD/kg 

Nu£q heat-transfer Nusselt number for drop of diameter D = D , 

hD 20 /k g 

n number of drops in given size range 

P c chamber pressure, lb/sq in. abs 

Q heat requirement, g-cal 

Re Reynolds number of single drop, D AVp /p 

§ gi 

Re 20 Reynolds number based on drop of diameter D = DgQ, ^20 AVpg/|ig 
Sc Schmidt number, Hg/bg v 

T temperature, °C 

AT difference between gas temperature and surface temperature of 
drop, T g - T g , °C 

T-p average film temperature, (Tg + T s )/2, °C 
t time, sec 

V velocity, ft /sec 

AV relative velocity of combustion gas with respect to drop, cm/sec 
V 20 velocity of a drop in spray having diameter D = ft /sec 
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w 

total weight flow, lb/sec 


X 

distance downstream from injector face, in. or ft 


A 

evaporation constant, 4kg At/p^Hp, cm^/sec 


P 

viscosity, poises 


P 

density, g/cc 


a 

molecular diameter, cm 


Subscripts : 


c 

completely evaporated 


d 

drop 


f 

film or final 


g 

combustion gas 


i 

initial 


l 

liquid 


m 

maximum 


n 2 

nitrogen gas 


p 

propellants 


S 

surface 


stoic 

stoichiometric 


T 

total 


t 

nozzle throat 


Y 

vapor or volume 


y 

condition for evaluating temperature 


20 

based on mean drop diameter, D^q 


30 

based on mean drop diameter, D 
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APPENDIX B 


EVALUATION OF d(ln N T )/d(ln d| 0 ) 

The diameter of a drop that is just completely evaporated D c in 
an incremental distance Ax may be calculated from the expression 

D c = 

where th e a verage Nujsselt number Nu is based on the average Reynolds 
number Re = D(V g - Vj)pg/|jg, as given in the relation 

Nu = 2 + 2.58X10 6 Re^~^ ^ 

In the definition of Re, D = D c /2, V z is the average drop velocity over 

the distance Ax, and V g is the average gas velocity over the same dis- 
tance. 

By means of the following trial-and-error method, the complete 
vaporization of small drops from stations x of 6 to 10 inches may be 
determined. If it is assumed that D c — 60 microns, then D — 30 mi- 
crons. Since at x = 8 inches V g = 268 feet per second and V^ = 145 
feet per second, AV = 123 feet per second. Therefore, Re = 18.9, 

= 2.91, D c = 49.6 microns. The assumption that D^ — 60 microns was 
too high and the calculation is repeated for D c - 40 microns, which 
results in D c = 46 microns. Finally, it is assumed that D c — 50 
microns, which yields D c =* 48 microns. Therefore, D c = 49 microns is 
taken as the final value. 

From the drop-size distribution data at x = 6 inches, a plot of 
D=D /D=D m 

S n / D n as a function of D was obtained (fig. 15). For this 
D=0 / D=0 

condition D = D c = 49 microns and the fractional number of drops in 
the spray that are completely vaporized is 0.10. Also, since 

D=D C 

E " 

ANgn = Niji — — — = 516(0.10) = 52 drops 
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Then at x = 10 inches 

(N t ) = (N t ) - £N t = 516 - 52 = 464 drops 

The same method of calculation is repeated to determine Np at 
x = 15 inches. The values of Nip are plotted against values of d|q, 
obtained at the same station x, as shown in figure 16. From this plot, 
it is evident that 

d(ln N t ) 

— = 0.45±0.02, 

a(m d| 0 ) 

which agrees within ±5 percent of the evaluation made in reference 1. 
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TABLE I. - DROP-SIZE-DISTRIBUTION DATA 





TABLE II. - FILM PROPERTIES EVALUATED AT AVERAGE FILM TEMPERATURE 
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Figure 1. - Diagram of experimental apparatus. 
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Figure 2. - Detail of propellant injector. 


Chamber pressure, P c , lb/sq in. abs 



(a) Fine spray. 



(b) Coarse spray. 

Figure 3. - Chamber-pressure fluctuations. 
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Figure 4. - Frequency spectrum of 
chamber -pres sure fluctuations. 
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Figure 5. - Photomicrograph of fine spray taken 10 inches downstream 
from injector face. X14. 


Drop velocity, V., ft/sec 


ro 
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Drop diameter, D, microns 


(a) Fine spray. 

Figure 6. - Drop-size and drop-velocity data. 
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Drop velocity, V=, ft/sec 
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(5) Coarse spray. 


Figure 6. - Concluded. Drop-size and drop-velocity data. 









Experimental Nusselt number, Nu. 


Fuel spray- 
Fine Coarse 


Evaluation temperature at - 


Experiment al 
Nusselt number, 


Experimental 
Reynolds number, 
Re 20 


+ T s )/2 
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Hr, 

g-cal/g 
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T f = (T g + T 
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T 

stoic 

T f ” ^stoic 
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*^stoic 

* ^stoic 


[1 

T s 

T e 

V 

O 

T s 

t e 


7 

Equation (7) 

I J 


+ T b )/2 



12 16 20 
(Re 2Q Sc gl/c 2 )°- 6 (k /k ) 0 - 5 


Figure 9. - Correlation of experimental Nusselt number with equation (7). 
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Figure 10. - Relation between velocity and diameter of 
mean drop size. 
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Figure 11, - Relation between total number of drops 
in spray and mean-diameter. 




Volume-mean drop diameter cubed, D (microns) 
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Figure 12. - Vaporization rate of volume-number mean drop diameter. 






Downstream distance, x, in. 


Figure 14. - Volume of fuel spray evaporated in 
rocket combustor. 





Drop diameter , D, microns 


Figure 15. - Relation between 


D=D 


B-H, 

E " 


D=0 / D=D 


and drop diameter 


for coarse spray; distance 
downstream from injector 
face, 6 inches. 




Total number of drops in spray 
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Area-mean drop diameter squared, 

d| 0 , (microns) ^ 

Figure 16, - Relation between total number of 
drops in spray and area-mean drop diameter 
squared. 
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